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Abstract. Invariant linearization criteria of square systems of second-order quadratically 
semi-linear ordinary differential equations (ODEs) that can be represented as geodesic equa- 
tions are extended to square systems of ODEs cubically nonlinear in the first derivatives. It 

O ■ is shown that there are two branches for the linearization problem via point transformations 

for an arbitrary system of second-order ODEs. One is when the system is at most cubic in the 

^ ■ first derivatives. We solve this branch of the linearization problem by point transformations 

j^ ■ in the case of a square sytem of two second-order ODEs. Necessary and sufficient conditions 

for linearization by means of point transformations are given in terms of coefficient functions 
of the system of two second-order ODEs cubically nonlinear in the first derivatives. A con- 
sequence of our geometric approach of projection is a re-derivation of Lie's conditions for a 
single second-order ODE and sheds light on more recent results on them. In particular, we 
show here how one can construct point transformations for reduction to the simplest linear 
equation by going to the higher space and just utilising the coefficients of the original ODE. 
We also obtain invariant criteria for the reduction of a linear square system to the simplest 
system. Moreover, these results contain the quadratic case as a special case. Examples are 
given to illustrate our results. 



1. Introduction 



A linearization problem involves the study of families of equations that are reducible via 
admissible transformations, which can be point, contact or more general, to linear equations. 
Lie [T] presented linearizability criteria, obtaining both algebraic and practical criteria, for a 
single second-order ODE to be point transformable to a linear equation via invertible changes 
of both the independent and dependent variables. 

Lie [1] proved that necessary and sufficient conditions for a second-order ODE, y" = E{x, y, y'), 
to be linearizable by means of invertible point transformations are that the ODE be at most 
cubic in the first derivative, viz. 

y" + Es{x,y)y'' + E2{x,y)y'^ + E,{x,y)y' + Eoix,y) = (1) 

and the coefficients Eq to E^ satisfy the over-determined integrable system 
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bx = —i^Eiy + -E2x + be — EqE^, 

by = E:,, - b^ + bE2 - E^Es + eEs, 
Ex = Eoy + e^ - eEi - bEo + E0E2, 

Cy = -E,y ~ -E2x - be + E0E3, (2) 

where b and e are auxiliary variables and the suffices x and y here and hereafter refer to 
partial derivatives. Since the classic work of Lie there has been continuing interest in this 
topic. We, inter alia, re-derive the Lie conditions ([2]) geometrically, by projections. 

Tresse [2] also studied the linearization problem for scalar second-order ODEs. He deduced 
two relative invariants of the equivalence group of point transformations, the vanishing of 
both of which gives necessary and sufficient conditions for linearization of equation ([T]). These 
conditions are equivalent to the Lie conditions (2) (see Mahomed and Leach [5]) and can be 
given as the compatibility of (2) as 

3(-Ei-E3)x — Eiyy + 2E2xy — 3{EoE3)y + E2Eiy — 2E2E2X — 'iEsxx — ^E^E^y = 0, 
3(-Eo-S3)x + 2Eixy — SEoyy " E2XX " E1E2X + 2EiEiy — 3{EoE2)y + S-Eo-^sx = 0. (3) 

Note that under the interchange of E-^ by ~Eq, E2 by —Ei and x by y, these conditions imply 
each other. Equations ([3]) provide practical criteria for linearization of equation ([1]) by point 
transformations. These conditions were also derived by the Cartan equivalence method (see 
Grissom et al ^ ) as well as recently using a geometric argument in Ibragimov and Magri [5] . 



The reader is also referred to the review of various approaches in Mahomed ^. Linearization 
via point and other than point transformations is of great interest and has been investigated 
in several works (see, e.g. [71[8l[9l[T0l[ni[T2l[l3l[Tl]). 

The algebraic criteria of linearization of systems of second-order ODEs by means of point 
transformations have been considered in Wafo and Mahomed [llj. Practical criteria for 
quadratic semi-linear systems of second-order ODEs have been looked at recently as well (see 
Mahomed and Qadir [IS]). In this paper our intention is to extend these results to cubically 
semi-linear square systems of second-order ODEs using geometric methods developed earlier 
(see Feroze et al [16]). As a by-product of our approach we re-derive the Lie conditions 
([2]). Moreover, we present practical criteria in terms of coefficients for cubically semi-linear 
systems of second-order ODEs to be linearizable by point transformations. As a consequence 
we provide practical criteria for the class of linear second-order system of two ODEs to be 
reducible to the simplest system. Notwithstanding, our results subsume the linearization 
criteria for the quadratic case. 

The outline of this paper is as follows. In the next section we present mathematical prelim- 
inaries. In section 3 we give an alternative method for obtaining the Lie conditions ([2]) as 
well as an alternative method for the construction of linearizing transformations for scalar 
second-order ODEs. Then in section 4 we derive practical criteria for linearization for a 
system of two second-order cubically semi-linear ODEs. Herein we state the relevant result 
for linear systems. Our theorem also contains the quadratically semi-linear equations as a 
corollary. In the next section we provide examples that amply illustrate our results. Finally, 
in section 6 we present a brief summary and conclusion. 



2. Preliminaries 

We first present some preliminaries. The system of geodesic equations is 

+ T)^x'x'' = 0, i,j,k = l,...,n, (4) 



x' 



where the dot refers to total differentiation with respect to the parameter s and r*j^. are the 
Christoffel symbols, which depend on x* and are given in terms of the metric tensor as 

r}fc = -g^"^{gjm,k + gkmj — gjk,m)- (5) 



The Christoffel symbols are symmetric in the lower pair of indices and have n'^{n + l)/2 
coefficients. The Riemann curvature tensor is 

r}i -pi -pi 1^ -pi -pm -pi -pm fn\ 

^jkl~^jl,k ^ jk,l ~^ ^ mk^ jl ^ ml^ jk^ \^) 

which is skew-symmetric in the lower last two indices and satisfies 

R)ki + R\ij + R\jk = 0- (7) 

A necessary and sufficient condition for a system of n second-order quadratically semi-linear 
ODEs for n dependent variables of the form (jl]) to be linearizable by point transformation 
and admit sl{n + 2, R ) symmetry algebra is that the Riemann tensor vanishes ([9l [IZ]), i.e. 

R)ki = 0. (8) 

Practical criteria and the construction of point transformations are given in [TS]. In particular, 
for a system of two geodesic equations (jl]), one has the linearization conditions (admittance 
of s/(4, R ) symmetry algebra) on the coefficients given by 

tty — b^ + be — cd = 0, 

by — Cx + (ac — 6^) + (bf — ce) = 0, 

dy — Cx — {ae — bd) — {df — e^) = 0, 

ib + f)x = ia + e)y, (9) 

where the Christoffel symbols are 

r\, = -a, r}2 = -b, r^2 = -c, rl, = -d, rl, = -e, rj^ = -/. (lo) 

Now equation (^ together with flTUl) on setting gu = p, gi2 = Q = fi'21 and (722 = 1^ yield 

Px = -2{ap + dq), 

qx = —bp — (a + e)q — dr, 

Tx = -2{bq + er), 

Py = -2{bp + eq), 

qy = -cp - (6 + f)q - er, 

Ty = -2{cq + fr). (11) 



The construction of the linearization point transformations are found as follows (see 
One invokes 
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where x = [x^, . . . , x"), u = [u^, . . . , u") with the requirement that gij{u) be the identity 
matrix. For the case of two variables, we need to solve the equations 

ul + vl= p, UxUy + v^Vy = q, ul + vl = r, (13) 

for which we have set (x^x^) = {x,y), (m\m^) = {u,v), gn = p, Qu = Q = fi'21 and (722 = r 
in (m. 



Following Aminova and Aminov [T7], we project the system down by one dimension and write 
the geodesic equations (jlj) as 

x'"' + Ak^x^'x^'x^' + B^x^'x^' + C^x^' + 0" = Q, a = 2,...,n, (14) 

where the prime now denotes differentiation with respect to the parameter x^ (in [1^ x^ is 
used as the parameter) and the coefficients in terms of the F^^s are 

A, = -FL, Bl = Tl - 25f,Fi)„ C^ = 2Ft, - 5^,T\,, D» = F?,, a, 6, c = 2, . . . , n, (15) 

where we have used the notation T(^a,b) = {Tab + Tba)/^. It is straightforward to deduce fin|) 
and (fT5|) . Indeed, insert 

.a_dx^.i _^ 

dx^ 
and its derivatives 

t// Ty tij I ^ ^^ • (X ^. .... It 

dx^ dx^ 

into system (j4]). These, after cancelation of x^ , directly yield (fT^ and (TT5|) . Note that in 
projecting down the Christoffel symbols there is degeneracy which results from the reduction 
of the range of the indices, so that T\^ and F^^^ appear in the same combinations in C^ and 
B^^, respectively. Consequently the set of coefficients A, B, C, D have n less elements than 
the coefficients F*;,. 



3. Re-derivation of the Lie conditions 



We invoke equations (TT^ and (TT5|) for n = 2. We also use (fTOl) in identifying the F*;.s with 
the coefficients a to / of the system of two geodesic equations which projects to (fT4l) . Thus 
we have (setting (x^,x^) = (x, y)) 

y" + Esix, y)y'^ + E^ix, y)y'^ + Ei(x, y)y' + Eo(x, y) = 0, (16) 



where 



E, 


= M2 = -r22 = c, 




E2 


~ -'-'22 ~ -'- 22 ~ ^-'^ 12 ~ 


-f + 2b. 


El 


_ n^ — opi r^ — 

— (^2 — -^-L 12 -^ 11 ~ 


-2e + a, 


Eo 


= D^ = Tl, = -d. 





(17) 

To re-derive the Lie conditions ([2]), we use the system of two geodesic equations (jl]) from which 
equation fITB]) arises projectively. Hence we utihze the conditions (jH]) which are conditions for 
a flat space. This requires that the coefficients a to / be in terms of the EiS. From (TT7j) we 
have 

a = Ei + 2e, 

c = Es, 

d = -Eo, 

f = 26-^2, (18) 

where we have chosen b and e as yet arbitrary. These are constrained by the relations iQ. 
We substitute ( TTS!) into iQ. Equations iQ then yield 

Ely + 2ey -b^ + be + E^Eq = 0, 
by - E3. + EiEs + eE3 + b^- bE2 = 0, 
Eoy + e^ + eEi + e^ - bEo + EqE^ = 0, 

36, - 3ey - E2, - Ely = 0. (19) 



The first and last equations of (TT9|) are easily seen to be equivalent to 
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bx = ^i^Eiy + -E2X — fe — -E0-E3, 



ey = ^E2x-'^Eiy-be-EoE3. (20) 



The second and third equations of (1191) as well as equations (l20l) . on replacing e by — e, are 
precisely the Lie conditions ([2]). Hence, we have provided an alternative derivation of the Lie 
conditions ([2]) by viewing the projection flT6|) in one higher space and looking at the fiat space 
requirement there. If we had projected the system of two geodesic equations to a single ODE 
of the form ( IT6l) by using x^ instead of x^, then by interchanging E^ by —Eq, E2 by —Ei and 
x^ = X by x^ = y, the coefficients ( TT71) imply the coefficients of the projected equation with 
independent variable x^. We state the following theorem. 



Theorem 1. A necessary and sufficient condition that the scalar second-order ODE flT6|) has 
s/(3, ]R ) symmetry algebra is that there is a corresponding system of two geodesic equations 
of the form (JH) from which it is projected that admits the sl{4, IR ) symmetry algebra. 



Furthermore, one can construct linearizing point transformations for flT6|) that satisfy ([3]) 
by resorting to the corresponding system of two geodesic equations from which flT6|) arises 
by projection. This is done by using the relations (IT3|) . This approach also results in the 
determination of at least one metric as a bonus. Notwithstanding, this method uses the 
coefficients of the equation which is linearizable and a transformation is then constructed via 
the relations (TT3|) . We consider two examples to illustrate this. 

1. On using flTSl) . the simple nonlinear equation 

y" + y'^-y' = (21) 

has corresponding a to f values, 

a = -1 + 2e, c = 1, d = OJ = 2b. 

These together with the choices 6 = and e = 1 satisfy the system (Q. With these values of 
a to / we obtain from (ITT]) particular solutions for p, q and r given by 

p = r = exp(2y — 2x), q = — exp(2|/ — 2x). 

Invoking (fT3l) . a linearizing point transformation to the simplest second-order ODE is 

?i = — =exp(-x + |/), t; = ^exp(-x-y), 

where u is the new independent variable. 

2. The familiar nonlinear ODE (see, e.g. [18]) 

y" + 'iyy' + y^ = Q 

has, upon using (ITS]) . 

a = 3y + 2e, c = 0, rf = -y^ , f = 2b. 

These and the choices b = 1/y and e = —y satisfy ([9]). A particular solution of (ITT]) is then 

p = 1 + X — 2xy^ + y^ , q = [l + X )y^ — xy^ , r = y^ {I + x ). 

A point transformation that linearizes the original ODE to the simplest second-order equa- 
tion, after solving (fT3]l . then is 

-X^ . 

2 !/■ 
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where u is taken as the new independent variable. This transformation was previously ob- 
tained in [IS] by mapping generators to canonical forms. As such we have presented another 
way of finding such transformations. 



4. Linearization conditions for square systems 



Driven by the success in obtaining the Lie conditions ([2]) by projection and then going back 
to the geodesic equations, we pursue similar conditions and practical criteria for linearization 
for a system of two second-order ODEs in a similar manner. Consequently, we study flT^ for 
linearization via point transformations by resorting to a system of three geodesic equations 
dl]). Before we do so, we need to first understand what is meant by linearization for systems 
of ODEs. A system of two second-order linear ODEs can possess 5, 6, 7, 8 or 15 point 
symmetries (see fi9[ 120]). The maximal symmetry algebra and hence sZ(4, R) is reached 
for the simplest system. Here we consider practical linearization criteria in terms of the 
coefficients for a system of two cubically semi-linear second-order ODEs of the form flT^ 
having s/(4, IR) symmetry algebra. The quadratically semi-linear case was treated in [15] . 
Also algebraic criteria for systems of second-order ODEs have been found in [11] . 



We once again invoke equations flMj) and flT5|) but now for n = 3. We therefore have 

x'" + A,,{x''f + 2A,,{x'yx'' + A,,x'\x''f + BUx^'f + 2Bl,x'' x'' 

+BUx^'f + Clx""' + Cfx^' + D^ ^ 0, 

x^" + A,^{x^'fx'' + 2A^,x^\x''f + A,,{x''f + B^x^'f + 2Bl,x^' x'' 

+BUx'''f + Clx^' + Clx^' + D^ = Q, (22) 

with coefficients 

A,, = -rL, Bl = Tl - 251TI),, CI = 2T1, - 5^,T\,, D^ = T^,, a,b,c = 2, 3. (23) 

Here three F^^ coefficients are lost. We select TI2, F^g ^^d F33 as arbitrary. We solve for the 
15 F^^s of ([23D in terms of the 15 coefficients Abe, B^^, Q, D^ as well as Tl^, Tj^ and F^g. 
We only write down the F^^s in which the arbitrary elements appear. They are 



11 - 


- ^J- ^2 '-^2 5 


:1 

13 - 


= 2*^^33- ^33) > 


<2 

22 - 


= ZL ^2 ' -022) 



T-2 _ 
■■-23 - 


= 2*^^33 + 2523-533), 


p3 _ 
J- 13 - 


- ■■- 12 + 2 3 2^' 


T-3 _ 
■■-23 " 


= i 12 + -D23, 



(24) 



The others can be read-off from equations ( l23l) . 



The ffat space requirement for the corresponding system of three geodesic equations (jlj) are 
now imposed by means of the vanishing of the Riemann tensor, viz. ([8]). They are (let 

{x^, x^,x^) = {x,y, z)) 

f-pi \ f-pi \ _[_ -pi -pm -pi -pm rj 

{'- j2)x {'- jlJy ~r i- ml'- j2 '-m2'-jl—^^ 

f-pi \ /-pi \ I -pi -pm -pi -pm r\ 

K'- j3)x K'- jl)z -r i- ml'- j3 '-m3'-jl~^^ 

U j3)y ~ i'- j2)z + J- m2'- j3 ~ '^ m3'- j2 ~ U' l^OJ 

which provide 27 conditions. Only 24 of them are linearly independent due to the identity 
(ITj). The reduction of these equations to explicit form is given in the Appendix. 

These are 24 conditions (147|) to (H9|) given in the Appendix that arise from the vanishing 
of the Riemann tensor as given in fl25|l . They are the Lie-type integrability conditions for 
the Fk. We find that there are 7 equations in P7|) to (H9l) which are independent of the 
r*;,. The other 17 contain first-order partial derivatives of the F*;.. Of these, F^2,y ^^^ ^i2,z 
appear once each, F33 ^ occurs three times and the rest twice each. Therefore, apart from the 
7 conditions which are independent of the FV and given solely in terms of the coefficients of 
the system, there arise a further 8 conditions on the coefficients upon equating the respective 
F*^. Hence, we end up with 15 conditions or constraint equations on the coefficients. Now 
the F*^ which appear once each do not result in linearly independent equations as can easily 
be checked by equating them with the corresponding F^ that were discarded. The resultant 
two equations that occur in this manner are linearly dependent. Thus the Ffg,?/ ^^^ ^u.z ^^^ 
spurious. It is thus opportune to state the following theorem. 

Theorem 2. A necessary and sufficient condition for the system of two cubically semi-linear 
ODEs 

y" + A,,y'^ + 2A2sy"z' + A,,y'z'' + Bl,y'' + 2Bl,y' z' + 53^3^'^ + Cly' + Cjz' + D' = 0, 
z" + A22y'^z' + 2A23y'z'^ + A332" + B^^y'^ + 2Bl^y'z' + Bl^z'^ + Cly' + C|z' + D^ = 0, 

(26) 

(where the prime denotes differentiation with respect to the independent variable x and the 
coefficients are in general functions of x, y, z) to be linearizable via point transformations to 
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the simplest system of two second-order ODEs is that its coefficients satisfy the following 
fifteen conditions on the coefficients functions of fl2BD. viz. 



\cl -Dl + \clCl + \clCl - D^Bl, - D'Bl, = 0, 

-D22X - 2 2y ~ ^22i>' + -(--2-022 + 2 3-°22 " 2 2-°22 " 2 23<-^2 " U, 

114 2 2 1 

-S23X - o-^22x + 2^*2?/ - o-^ ^23 - 2-B22C3 + 0-^23^2 " ^^Sy = 0, 



Cl - Dl + tCld + tCld - Bl,D^ - Bl,D' = 0, 



r2 r"^ n2 /I 1/^2753 1 r2 ^2 r2 /^3 I 1 r2 ^2 _ n 

^33x - 2 3^ - />> /i33 + -03^^33 - -i^23'-'3 " 3 33'-^3 + 2^'ii^2 " U, 

—^23?/ + ^222 — ^22-^23 ~ ^23-^23 + ^23-^22 + ^33-822 = 0, 

-^33?/ + ^232 - ^22-833 - ^23-833 + ^23-823 + ^33-823 = 0, 

5 11 1 

-^23a' + 2^23^2 + 2^33^2 " 0-^232 + -^33-^22 + 7^*3^23 - -B23-B23 

2 12 1 

-3-^231/ + 3^332^ + 3-SI22 - -C'1^22 = 0, 

1 1 

-^33x + 2*^2^33 + 2^33^3 " -^33^ + ^23^ - 532^33 + 523^23 - ^23^33 + 533^23 = 0, 



q^22x + o^2y - 2^2-^33 + ^ A22 - -D A23 - -C^B22 + g-^23'-^2 



_j_ r3 /^3 1^ r^3 r)3 /^2 7^3 ri 

+-°23x ~ 2 22 + 2 3-°23 ~ 2 2-^23 " U, 

-^22x + -C2A22 - -B22-B33 + -B23J/ ~ -^222 + B22B23 + -823-^23 + 2^3^22 " -^23-^22 = 0; 

Dl + i?2^2^' + D'Bl, - D^Bl, + ^Cl - ^Cl - Dl + ^ClCl - \clCl - BI,D' = 0, 

-2^23^ + 3-833^ + -A23C2 + 3^23^3 + 3^3^22 - 3^232 " 3*^2^33 + 2522^33 

2 2 

-2523^23 - 3^23y + 3^222 = 0, 

r2 j_ r''^ 9 n2 /I I r'2 p3 I 1 ^2 r2 I 1 ^3 r2 ^ p2 ^2 p2 ^3 

^23x + 2 3y ~ ^3 2 3 23 + 2^3-^22 + 2^3-^23 ^ 2 23'-^2 " -D33<-^2 

-C^, - ^3^33 = 0, 

r2 _l r3 _i_ /^2 /^2r3 I /^2 r2 , r2 ^3 r2 ,^2 ^ ^3 

~-D23x + -D33X + "-"Sj/ ~ '-^3-'^23 + '-'3-'^22 + -'^23'-^3 ~ -'^23'-'2 ~ 9 32 

~Cl - 2D' A,, = 0. (27) 

Proof. The proof follows from the preceding discussions. For if the system of two equations 
f l26|) are linearizable by point transformation to the simplest system, then its coefficients 
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can be written in terms of F*^ as in equations fl23|) which in turn gives rise to the Lie-type 
integrabihty conditions on the P-^ and hence fl7r|) . Conversely, if the coefficients of the 
system of equations (126!) satisfy the fifteen constraint conditions on the coefficients given 
by the relations (j27j) which is a consequence of the Lie-type conditions P7|) to (149|) . then 
the coefficients of the system (126!) can be written in terms of the F';, and the corresponding 
geodesic equations in three-space is linearizable as well as the projected equations (12BI) . 

Corollary 1. The system of two quadratically semi-linear ODEs 

y" + Bl^y'^ + 2Bl^y'z' + 53^3/^ = 0, 

z" + Bl^y'^ + 252^3^'^' + Bl^z'^ = 0, (28) 

where the B'^^s are functions of y and z and the dot denotes total derivative with respect 
to X, is linearizable by point transformations to the simplest system of two equations if and 
only if the ^^^.s satisfy the four conditions on the coefficients given by 



d3 d3 I d3 d3 I d3 752 , 753 753 753 752 
~ ^22^33 + ^2'iy " ^222 + ^22^23 + ^23^23 ~ ^23^22 " 


= 0, 


753 753 I 9 753 752 9 753 752 752 _, 752 

o-°33j/ o-'-'23z "t ^-'^22-°33 ^-D23-°23 n 2'iy ^ o 'i2z " 


= 


1 2 12 

R3 I r2 753 R2 753 ^752 I -^ 753 1 ^ 752 

~i;^21.z "T -D33-D22 ~ -D23-D23 ~ o-D23y + 0-D33J/ + 0-D22Z " 


= 0, 


-^33y + ^23z ^ ^22^33 + ^23^23 " ^23^33 + ^33^23 = 


= 0. 


le sets B22 = —0-, -BI3 = —b, B^^ = — c, B22 = —d, -Bf3 = 


: — e and B^^ 



(29) 



-/, 

one gets precisely the conditions ([9]). Hence Theorem 2 naturally contains the linearizability 
criteria for the quadratic case. 

Corollary 2. The system of two linear (in the first derivatives) ODEs 

y" + Ciy' + Clz' + D' = 0, 

z" + Cly' + Clz' + D^ = 0, (30) 

where the prime refers to differentiation with respect to x and the C^s are independent of 
y and z, is linearizable by point transformations to the simplest system of two equations if 
and only if the C^s and D°s satisfy the three conditions on the coefficients, viz. 





^-ct^^-clcl^^-clcl -- 

2 23S 1^321^22 


= Dl 






2^3x + ^^3^3 + ']p-i^2 = 


- Dl, 




2^3. 


"" 2^2z + ^^'3(73 - -6*2^*2 = 


- D\- 


-Dl 



(31) 
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We have provided practical criteria, necessary and sufficient conditions, for equations of the 
form fl2B]) to be hnearizable via point transformations to the simplest system. The question 
naturally arises if there are more general equations than (l26l) that can be hnearizable to the 
simplest system. Indeed, there are more general systems of two second-order ODEs which 
can be linearized. 

The most general system of n — 1 second-order ODEs hnearizable is given by 

Jix^" + G^x'^'x^" + /^%ix^' x^' x'' + K%x''x^' + n)x^' + E* = 0, 2 = 2, . . . , n, (32) 

where the prime refers to total differentiation with respect to x^, the coefficient functions are 
dependent upon x^, . . . , x", and are given by 

J'i x^l \ri x^l x^i 

j — -^ 1^ j — ^ j^ li 

kj — ^,k^,j ^,j^,k^ 

jkl - ^,l^,jk - ^,jk^p 

E' = X\X],, - X\,X],, z, J, A;, / = 2, . . . , n (33) 

in which 

X^ =X^(x\...,a;"), X' = X\x^ ....x""), i = 2,...,n (34) 

are invertible transformations. It is certainly not difficult to obtain (132|) . This is done by the 
substitution of (1M|1 into the free particle system 

X^" = 0, ^ = 2,...,n; ' = ^. (35) 



This after routine calculations yields ( l32l) with the coefficients satisfying ( l33l) . Equation ( l32l) 
is the most general system of n — 1 equations point transformable to the simplest system 
( !35|) . Equation (132|) has n{n — l){Tn? + 6n — l)/6 coefficients. 

Equation fl32|) can be written in normal form in terms of at most cubic ffist order derivatives 



as 



x'" + A)^ix^' x''' x'' + Bi^x^'x^' + Cix^' + D' = 0, ^, J, /c, / = 2, . . . , n, (36) 



provided 

'-^klm — ■-'j^klm ^ ^mj-^kli 
^kl = J]Bkl + GljCf^, 

E' = Jp\ 



^pj^klm — 0- (37) 
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The relations fl371) can be obtained by solving for the second derivative in terms of the first 
order derivatives and inserting these into equation fl5^ . The last equation of fl57|) tells us 
that not all the A^.j^i^ coefficients are independent. As a matter of fact if we replace these 
by Aki in (!36|) . then it turns out that this relation in (!37j) will now be identically satisfied. 
What transpires is that the quartic term disappears automatically due to Gp^ being skew 
symmetric in the lower indices and x^'x^ appearing symmetrically. One also needs then to 
adjust the relation (lS7k) in the latter case by 

^\lm = Jk^lm + Gl^jBf^i. (38) 

The remaining equations of fl37p are the same. 

There are two branches of the linearization problem by point transformations for a system of 
n — 1 second-order ODEs. One is the general form (132!) owing to the arbirariness of the A* ^^.^ 
coefficients. The other is the form (fT4l) in which the cubic coefficients are fewer in number. 
In the case of two second-order ODEs, equations (I26p . we have obtained explicit linearization 
criteria as encapsulated in Theorem 2 and their corollaries. 



In the general equation (132!) there are (n — l)n(n^ + Qn — l)/6 coefficients while for ([T] 
there are {n — l)n{n + 2)/2 independent coefficients. It would be of interest to find practical 
criteria for the reduction of equation (152]) to the simplest system via point transformations 
for n = 3. Of course it is of great interest to do this for the general system ( I32p for n > 4. 

If one has a system of the form (!32l) with known coefficients which is reducible to the free par- 
ticle system (l35l) by point transformations, then one can utilise (l33l) to construct a linearizing 
point transformation. Also, we can obtain linearizing point transformations for system (TT^ . 
if it is linearizable to the simplest system (15^ . by invoking (155|) together with (ITTI) . 

In particular, one can find linearizing point transformations for the system (!26|) in a similar 
manner by solving the system ([3j 



Instead of using the system (|33|) in order to construct a linearizing point transformation there 
are other ways as pointed out earlier. One is to go to the higher space, once one has the 
coefficients at hand, and use (IT^ for which (7jj(u) must be the identity matrix and where 
we may set u^ to be the independent variable. Yet a third approach is that of mapping 
symmetry generators of the linearizable system, if known, to the free particle generators. 
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5. Examples 



We present examples to illustrate our results. We have y and z as the dependent variables. 
Also the ' below denotes differentiation with respect to x. Moreover, we have included one 
example that does not satisfy our linearization criteria but belongs to the more general class 
which is linearizable. 



1. Consider the anisotropic oscillator system 

y" + ui{x)y = Q, 

z" + uj2{x)z = Q, (39) 



The coefficients of system fl39|) satisfy the conditions fl3T|) provided 001 = 002- Hence in order 
for the system fl39l) to be reducible to the free particle system one must have isotropy. 

2. The simple linear system 

y" + z = 0, 

z" + z = 0, (40) 

do not satisfy the conditions flHT]) . Thus this system is not transformable pointwise to the 
free particle system. This system does not have a Lagrangian formulation as well [2T]. 

3. For the quadratic system 

y -y +y = o, 

z" -z' + z'^ = 0, (41) 

all conditions (l29l) are satisfied. Therefore the system (HT!) is reducible to the simplest system. 
A point transformation that does the job is 

u = expx, V = expy, w = exp^;, (42) 

where u is the independent variable. This can be constructed by going to the higher space 
as we have illustrated for the scalar ODEs in section 3. 

4. Consider the cubically semi-linear system 

y" + -y' + y" + i- + -2)y" = o, 
X y y 

z" + -z' + z'^ + 2y'z' + (- + ^)y'^z' = 0, (43) 

X y y 
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For the system fj43l) all the conditions fl27|) hold. A linearizing point transformation to the 
simplest system is 

u = Inxy, V = expy, w = exp{y + z), (44) 

in which u is the independent variable. 
5. Finally the system 

A 2 12 , A 2 / / I o 2 /3 I o /2 / I o 2 / /2 , o 2 / // 2 2 // , o 2 / // 

Ayz y + Ay zy z + 2xz y + Sxyzy z + 2xy y z + 2xy zy z = y z y + 2xy zz y , 

II X I II , I II 2 i2 II I I II 1/ I , i\/ I , I , r\ I I < li\ /'Ac:\ 

y + xzy y + xyz y — xz y y — xyzy zy = y [zy + yz )[zy + yz + 2xy z + xyz j, (45) 

is not of the form given in Theorem 2. It is of the form given in (l32l) and is linearizable by 
means of the point transformation 

u = xexTp{yz), v = xy'^z'^, w = y, (46) 

where m is the independent variable. 



6. Concluding remarks 

Aminova and Aminov [T7] had provided a procedure of projecting down 1 dimension from 
a system of n geodesic equations to n — 1 cubically semi-linear ODEs. Separately, we had 
provided [15] linearizability criteria for a square quadratically semi-linear system. These were 
used together to derive linearizability criteria for a single cubically semi-linear equation by 
projecting down from a system of 2 quadratically semi-linear equations. This provided an al- 
ternate method to prove Lie's general result for linearizability of a single non-linear equation. 
It led naturally to an extension of the linearization criteria via point transformations from 
a scalar second-order ODE as obtained by Lie [1] to a system of two cubically semi-linear 
ODEs of the form (!26l) . These provided necessary and sufficient conditions for reduction to 
the simplest system and hence s/(4, R) symmetry algebra for equations of the form (!26l) . 
Moreover, Theorem 2 provides criteria for the reduction of linear systems of two equations 
to the free particle system. 

Lie had demonstrated [1] that only cubically semi-linear scalar equations of order two are 
linearizable in general. As such, it could have been hoped that the projection procedure will 
provide the complete solution of the linearizability problem for the system of 2 non-linear 
ODEs. That hope is doomed from the start as there are 5 classes of systems of 2 cubically non- 
linear equations that are linearizable by point transformations, having different symmetry 
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algebras. Moreover, the maximum symmetry algebra class of such systems of two equations 
is one branch of the linearization problem via point transformations as the general class is 
represented by (l32l) . Why do we get a unique class in the former case and 5 in the latter? 
Furthermore, how many distinct classes should there be for a system of n cubically semi-linear 
ODEs? 

We start by noting that the projection procedure and linearizability can be equally well 
adopted for an arbitrary system of n quadratically semi-linear second order ODEs reduced to 
n — 1 cubically semi-linear second order ODEs. There are two branches for the linearization 
problem for systems admitting the maximal algebra for n > 3. There is enormous computa- 
tional complications that arise. As such, one would need an algebraic computational code to 
deal with larger systems. A code has, indeed, been prepared to construct the metric coeffi- 
cients given the Christoffel symbols [22] . That can be extended to deal with the linearization 
of larger systems. Now observe that in projecting down from the system of n dependent 
variables to n — 1 variables, the Christoffel symbols are reduced from n^(ra + l)/2 by n, to 
give {n — l)n{n + 2)/2 independent coefficients. Since we now have n — 1 equations, each with 
its own cubic function, there are (n — l)n/2 cubic coefficients for the reduced system. If the 
number of coefficients left over after losing n equals the number of coefficients of the reduced 
system, we can determine one set of coefficients in terms of the other. The two expressions 
are obviously equal for n = 2 and the former is greater than the latter for n > 2. As such, 
the coefficients of the cubic system can be determined uniquely in terms of the quadratic 
system for n = 2, i.e. for a scalar cubically semi-linear system. For larger systems there will 
be infinitely many ways to write the former in terms of the latter. Hence there is a unique 
solution to the linearizability problem only for the scalar cubically semi-linear equation and 
many solutions for systems of cubically semi-linear systems! 

The second question remains and has, in fact, been compounded. It is known that there 
are 5 and not infinitely many distinct classes. Why? The point is that all distinct ways of 
writing the cubic system coefficients in terms of the quadratic system coefficients will not 
give independent criteria as there will be transformations permissible from one definition to 
another. The point is to determine those that are distinct. Another way of looking at what 
we have done is to note that we have asked that the original system correspond to a system 
of geodesic equations in flat space. Then the projection gives the reduced system, which 
must also be of geodesies in an (n — l)-dimensional flat space. Even if the original geodesies 
were curved, the projected geodesies could correspond to straight lines. For example, if the 
original space was a sphere and one projects along the plane containing the geodesic to a 
plane perpendicular to it, the resulting projected curve would be a straight line. 
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The minimal dimension for a system of n second-order ODEs to be linearizable by point 
transformation is 2n+ 1. The maximum dimension of the symmetry algebra is (n+ l)(n + 3) 
which corresponds to s/(?2 + 2,1R). The other submaximal symmetry alebras besides that 
of dimension 2n + 1 range from 2n + 2 to (n + 2)^/2 for n even and [(n + 2)^ + l]/2 for n 
odd. Thus for n = 2 we have the mimimum dimension to be 5 and other submaximal algebra 
dimensions are 6, 7 and 8. The maximum dimension for n = 2 is 15. For n = 3 the minimum 
dimension is 7 and the next to maximum is 13. The maximum is 24. Thus for this case there 
are 8 classes. Generally, for n = 2m, the number of classes is 2r?7,^ + 3 and for n = 2m — 1 it 
is 2m^ — 2m + 4. 

It would be important to find ways of providing the linearizability criteria for the cases of 
the other symmetry algebras. 



Appendix 



We take j = 1 in the third set of fl25|) as the 3 dependent equations and discard them. The 
invocation of the first set of 9 equations of fl25l) gives 

Icl -Dl + ^ClCl + )clCl - D'Bl, - D^Bl, = 0, 

-D22Z - 2 2s/ ~ ^22^^ + 2'-'2-D22 + 2^3-^22 " 2 2-°22 ^ 2 23<-^2 " U, 

111 
^I2,y = —^22x — ^22ri2 + ^"2^22 + ^ V1^22 + ^121^12 + 2-^22^33 — 2-^22-^33 + 2^2^23, 

12 11 

p2 ___d2 I _r'2 _|_ pi r2 _|_ _/^3p3 __/^3p3 
■■- 12,y ~ o-°22x "T o^2y "r ■■- 12-'- 12 "r a^2'- 33 ^'-^2-'333 

2 11 

+D A22 + -D A23 - -C^B22 + -B23C2, 



9 1 11 

1 _ __r2 I r<2 I -pi -r2 I -^^3-r3 -^^^3^3 I n2 

2,x — o 22; 

1„. . 1 



^12,0; — ~o-^22z + o^2j/ + ri2ri2 + ^(^2^33 - 2^2-833 + D A22 



+-D^A23 - -CiB22 + --B23CI 

1 14 2.2 1 . 

-^232: - 2-S22a' + 2^*2^ " o"^ ^23 " 0-822^3 + 0-^23^2 " 2^3y = 0' 
^33,2/ = -2yl23a; + -633^^ - 2r^2^23 + ^23^2 + 2ri2-B23 + ri2r33 - r^2-^33 

+ ^33-^23 ~ -^33-^23 + ^22^3 + ^23^3 

p3 f) 752 I T53 I /^2 I r) 7~)2 A /^2 T53 I /^2pl p2 753 1 r^2 lYl , 02 /~i3 

J- 33,a- - -^-D23x' + -D33:r + '^3y + ^'-' ^23 " <--3-D23 + '-'3^^ 12 " ^^ 12-^33 + <-^3-D22 + -D23'-'3 
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r2 ^2 /^S r3 I d3 ,^2 I -*- ^3p3 ^^2^3 i pS pS {A7\ 

~^2'i'^2 ~ 2^3-033 + 2 33'-"2 + 2^ ^^ ^ 2 33 ~^ 33^- 12- l^' J 

The second set of 9 equations of fl25|) yields 



2<^2^ + -D ^23 - 2*^2^23 + |C'2^33 + |* 



pi oS I /^3 I r-)3 /I r^3 752 , /^3 753 , /^3t-i3 

J- 12,a; - -^22,x + 7,'-'2z + ^ ^23 - ^'-'2-023 + 7'-'2-D33 + 7<-^2J- 33 



/^3 r3 i^ /^2 753 j^ p2 pi 
~2 3-°23 + 2 2-°23 + ^ 12-'- 12) 

^12,2 = ~M?,x - ^23ri2 + ^23^2 + ^12-^23 + -r33S23 

1111 
p3 pi 753 753 753 pi 1/1 /^3 

2^ 33^^ 12 ~ 2 33-^23 ~ 2^^ ^^ 2 ^^ ^ ' 



^C7|. + ^C'L + Dl + 11^3^3 + \d^BI, - ^c _ , ^ . 



r?2,. = —A + 0C2I + ^' + ^rL^' + o^'^l3 - -Acl - -7CICI 



p2 _ r2 ,9/1 n2 r'2 p3 I 1 r2 ^3 I ^2 _i_-'-/^2pl , ''^ r^3p3 

J- 12,^ - --D23X + ^^23-'^ - 2^3-^23 + '^^Zz'^2 + '-^2z + 2^3^- 12 + ^'-'3^- 33 

r'2p3 I -'-r3 r2 03/^3 1 1^2p3 d3 p2 _, /i n3 

-^<-^2J- 33 + 2^ 33J- 12 - ^^33<-^3 + ^<-^2^33 ~ ^^^ ^^ ' 



^c^L - ^.^ + \clcl + \ 



Cl - Dl + -ClCl + -ClCl - Bl,D^ - BID' = 0, 



r2 1^2 p,2 A j_ n'iTyS -'- r2 ^2 r2 r^3 I -'- r2 ^2 _ n 

^33x - 2^3z - ^ ^33 + 2'-^3^33 " 3 23'-^3 " 3 33'-^3 + 2^33'-^2 " U, 

^33,^ = -S33X - 4523^. + 6A23D - 2C3523 + 2533C2 + 26*2^ + C3ri2 + 2C'3r33 



_^2p3 _i_ r3 r2 _ _ r3 r'3 I ^ 
'2^2-"- 33 "T -■- 33-"- 12 2 33'-^3 "T 2 



/^2p3 j^ p3 p2 733 /^3 j^ /^2 733 733 p2 1^ o /i r^3 

'-'2-'- 33 + -*- 33-*- 12 ~ o-°33'-"3 + o'-'2-D33 ~ -033-*- 12 + ^^33-'^ , 



p3 /^3 I /^2 752 I r) A 7-)2 I /^2pl I ,^3p3 /^2p3 

J- 33,x — 2*-'32 + 2^2z ~ -D23X + ^^23-t^ + ^3^ 12 + 3^ 33 " O ^^ 

I p2 p3 /^3 d3 I /^2 753 d3 p2 , o 4 7-)3 

+ i 12^ 33 - 2 3^33 + 2*^2^33 " ^33^ 12 + 2/i33l>' , 

1 1 

^33,2 = -2yl33a; + i?33^ - 2^331^2 + C2A33 + 2r^2-S33 + 2^33^33 ~ 9 "^33-^33 

+A23Ci + A33Cl (48; 

The last set of the 9 equations of (!25|) result in 6 independent conditions 

pi d3 I 753 I /^3 A 753 752 I 753 753 i 753 p3 

J- 12,y — ~^23y + -^222 + 9 '-'2 ^23 " -D22-D23 + 9 22-°33 + 9 22 J- 33 

"-^23-^23 + ri2ri2 — -C^A22 + -C2A22 — v422r;^2 + -^23-^22 + -^22^12; 

^12,2 = 0-^232 + gC'2v433 - -B33-B22 - g^3^23 + 2^*2^23 - ^23ri2 + -B23-B23 ~ 2 "^23-^33 



1 112 12 1 

+ 2^2?.^ 33 + -^23^- 12 ^ 2^'^ ^^ 2 ^^ ^^ Q 232/ ^ o-°33j/ ~ o-°22z + o'-'3^22, 

4 2 2 2 

^33,?/ = 0-^232 + 0^2^33 " '^B^^B22 - ^^^3^23 + ^^2^23 " 2^231^2 + 2-023-^23 

2 12 1 

"-^23-^33 + -^23^33 + 2-023^12 ~ -^33^12 + ri2r33 + ^-B23j^ — o-^33j/ ~ o-^222 + o^3^22, 

p3 o r)2 r) 752 I r)3 o /i p2 i r) 752 752 1 r) T52 pi I p3 p3 

J- 33,2 — ^-°33y ^ ^-D232 + -^332 " ^^33J^ 12 + ^-D22-D33 + ^-^33^^ 12 + g 33^ 33 

+C3A23 - 2i?23i?23 + 2-823-^33 ~ 2'^33-^33 ~ 2^33^23) 

—^23?/ + ^222 — ^22-^23 ~ ^23-^23 + ^23-^22 + ^33-^22 = 0? 

-A33j^ + A232 - ^22^13 - ^23^13 + ^2352^3 + ^335^3 = 0, (49) 
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